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Abstract 

We report a measurement of the charge asymmetry for same-sign dileptons in B^-B^ mixing. 
The data were collected with the Belle detector at KEKB. Using a data sample of 78 fb~^ recorded 
at the T(4S') resonance and 9 fb~^ recorded at an energy 60 MeV below the resonance, we measure 
A,i = (-0.13 ± 0.60(stat) ± 0.56(sys))%. 

PACS numbers: 13.65.+i, 13.25.Gv, 14.40. Gx 
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I. INTRODUCTION 



The standard model allows CP- violation in B^-B^ mixing 0]. In particular, there is a 
possible difference between the B^ — > B^ and B^ — B^ transition rates, which can manifest 
itself as a charge asymmetry in the same-sign dilepton events in T(4S') decays when prompt 
leptons from semileptonic decays of neutral B mesons are selected. With the assumption of 
CPT invariance in the mixing, the flavor and mass eigenstates of the neutral B mesons are 
related by 

|fi^)=p|50) + g|50), 

\BL)=p\B')-qm. (1) 
The time-dependent decay rate for same-sign dileptons is given by 

(2) 



\ (4:3)^1+ i+ {At) = -z-^e 

OTbO 



4 

Atl/r^O 



p 


2 






cosh (^^At^ 


Q 







for the sample. For the sample, p/q is replaced by q/p. Here Arriii and AF are 
the differences in mass and decay width between the two mass eigenstates, r^o is the average 
lifetime of the two mass eigenstates. At is the proper time difference between the two B 
meson decays. In this analysis only the absolute value of At is used. It is assumed that the 
semileptonic decay of the neutral B meson is flavour specific and CP conserving, so that 
Ai = Ai. If CP is not conserved in mixing, the condition |p/g| = 1 is no longer required 
and the decay rates for i^i^ and £~i~ samples can differ. As can be seen in Eq. |21 the 
At dependence is the same for the i^i'^ and £~i~ samples, and therefore the CP- violation 
shows up as a At-independent charge asymmetry, defined as 

^ ^ ^T{4S)^e+e+ - ^r(4S)^e-e- _ 1 - Iq/pl'^ ^ 4Pe(eg) 

Here corresponds to the ex describing CP- violation in the neutral K meson system. 
Standard Model calculations give the size of this asymmetry to be of the order of 10~^ 
A significantly larger value would therefore be an indication of new physics. 

Experimentally, measurement of same-sign dilepton events that originate from B^B^ and 
B^B^ initial states requires careful charge- dependent corrections, which must be done in 
several steps. First, the contribution from continuum e'^e~ — > qq (where q = u,d, s or c) to 
same-sign dilepton events must be subtracted using off-resonance data. Second, all detected 
lepton tracks must be corrected for charge asymmetries in the efficiencies for track finding 
and lepton identification, and for the probabilities of misidentifying hadrons as leptons. 
After these corrections, the remaining same-sign dilepton events still contain backgrounds 
from B^B^ and B^B~ events. The last step of this analysis is to separate the signal events 
from these background events using their different behavior in the At distributions. 



II. BELLE DETECTOR 

The data were collected with the Belle detector j3| at the KEKB asymmetric e^e~ col- 
hder Q. 
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The Belle detector is a large-solid-angle magnetic spectrometer that consists of a three- 
layer silicon vertex detector (SVD), a 50-layer central drift chamber (CDC) for tracking, a 
mosaic of aerogel threshold Cherenkov counters, time-of- flight scintillation counters (TOF), 
and an array of CsI(T£) crystals for electromagnetic calorimetry (ECL) located inside of 
a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron fiux-return 
located outside of the coil is instrumented to detect mesons and to identify muons 
(KLM). The integrated luminosity of the data sample is 78 fb~^ at the T(4S') resonance 
( "on-resonance" ) and 9 fb~^ at 60 MeV below from the T(4S') resonance ( "off- resonance" ) . 

A. Track finding Efficiency 

The track finding efficiency is determined by analyzing a sample where simulated single 
electron or muon tracks are overlaid on hadronic events taken from experimental data. 
Lepton tracks are generated to cover the region of 1.2 GeV/c < p* < 2.3 GeV/c and 30° < 
^lab < 135°, where p* is the lepton momentum in the e^e~ center-of-mass (cm.) and 6'iab is 
the angle of lepton track with respect to the 2;-axis in the laboratory frame. The z-axis passes 
through the nominal interaction point, and is anti-parallel to the positron beam direction. 
Figure n shows track finding efficiencies for positive and negative tracks separately and their 
fractional differences as a function of p* for electron and muon tracks. Events in all 6'iab. 
regions are combined in these plots. The charge dependence of the track finding efficiency 
for both electrons and muons is less than 1.0%. 

B. Lepton Identification 

The most important contribution to the electron identification comes from examination 
of the ratio of the ECL cluster energy to the track momentum measured in the CDC. This 
information is combined with the shower measurement in the ECL, the specific ionization 
measurements {dE/dx) in the CDC, and the ACC light yield, to form an electron likelihood 

The two-photon process e+e^ {e'^e~)e~^e~ is used to estimate the electron identifi- 
cation efficiency. For this data sample, events are required to have: i) two tracks with 
particle ID information inconsistent with a muon hypothesis, laboratory momenta greater 
than 0.5 GeV/c and transverse momenta greater than 0.25 GeV/c; ii) at least one ECL 
cluster with energy greater than 20 MeV. The two tracks are required to have: i) an acolin- 
earity angle whose cosine is greater than —0.997; ii) a transverse-momentum sum less than 
0.2 GeV/c; and a longitudinal momentum sum of less than 2.5 GeV/c in the cm. frame; 
iii) an invariant mass less than 5 GeV/c^. In addition the sum of the ECL cluster energies 
must be between 0.6 GeV/c and 6.0 GeV/c. The electron identification efficiency is obtained 
by taking the ratio of the number of tracks selected with the above requirements with and 
without additional electron identification requirements. 

For muon identification, CDC tracks are extrapolated to the KLM and the measured 
range and transverse deviation in the KLM is compared with the expected values to form a 
muon likelihood Q. 

The muon identification efficiency is determined by analyzing a data sample where simu- 
lated single-muon tracks are overlaid on the hadronic events taken from experimental data. 

Figure El shows the charge-dependent lepton identification efficiencies, where electron 
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FIG. 1: Track finding efficiencies as a function of cm. momentum for positron tracks e+ (a), 
electron tracks e_ (b), and charge dependence defined as (£+ — £-)/{£^ + £_). Corresponding 
plots for muon tracks are shown in (d), (e), and (f) 

tracks arc required to satisfy Cf. > 0.8 and the muon tracks are required to satisfy > 0.9 
and the reduced x^, of the transverse deviation in the KLM is required to be less than 3.5. 
The charge dependence of both electron and muon identification efficiencies are less than 
1%. 

C. Hadron identification 

Comparison of the ACC light yield to the track momentum, the time of flight measure- 
ment, and dE/dx measurements in the CDC are combined to provide hadron likelihoods, 
vCtt for pions, Ck for kaons, and Cp for protons. 
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FIG. 2: Identification efficiencies as a function of cm. momentum for positrons e+ (a), electrons 
e_ (b), and charge dependence defined as (e-|- — e-)/(£^ + £_)(c). Corresponding plots for muons 
are shown in (d), (e), and (f). 



D. Fake lepton 

The hadron fake rate, which is defined as the probability that a hadron track is mis- 
identified as a lepton, is determined from a sample of Kg — > 7r"'"7r~ for pions, K^K~ 
for kaons, and A — > pn~ (A — > pn^) for protons. These decays are selected from a hadronic 
event sample, which will be described later. To select these track pair combinations, the 
closest approach with respect to the run-dependent interaction point and the position of 
decay vertex are used. The z position distance of two tracks and the deflection angle 
(except for — > KK) at the decay vertex are also used. For each decay, the invariant mass 
of the two tracks is calculated after imposing the hadron identification requirement on the 
negative (positive) charged track. The signal yields are obtained by fitting the resulting mass 
distributions to sums of double Gaussian signal terms and smooth background functions in 
two ways: once without imposing any particle identification requirement and again after 
imposing the lepton identification requirement on the positive (negative) charged track. 
The ratios of the two signal yields give the fake rates for the positive (negative) charged 
tracks. The following cuts are placed on the likehhood ratios: Ct^/{Ct^-\- Ck) > 0.8 for pions 
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in Ks — > 7r"'"7r~, CkH^^k + -Ctt) > 0.8 for kaons in — > K^K^ , and Ct,/{Ct, + £p) > 0.8 
for pions in A — pTT~ . The rate of pions faking electrons is at most 0.1% for both charges 
and shows no significant charge dependence. The rate of kaons faking electrons decreases 
rapidly as piab becomes larger and is less than 0.2% for piab > 1.4 GeV/c with no significant 
charge dependence. While the rate of protons faking electrons is nearly zero, the rate for 
anti-protons faking electrons is as large as 4% due to the large anti-proton annihilation cross 
section in the ECL. The rate of pions faking muons is about 1% for piab > 1-5 GeV/c and 
shows no significant charge dependence. The rate of kaons faking muons is 1% to 2% and 
that for is about 50% larger than K' due to the larger kaon-nucleon cross section for 
the K~ . The rate of protons faking muons is less than 0.4% and shows no clear charge 
dependence. 
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FIG. 3: Rates of pions ((a) and (d)), kaons ((b) and (e)), and protons ((c) and (f)) faking electrons 

and muons vs laboratory momentum. Filled circles are for positive tracks and open circles are for 
negative tracks. The increase in the rate of kaons faking electrons at low momentum clearly visible 
in (b) is due to the overlap of the electron and kaon energy loss bands. 
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III. EVENT SELECTION 



A. Hadronic Event Selection 

Hadronic events are required to have at least five tracks, an event vertex with radial and 
z coordinates within 1.5 cm and 3.5 cm respectively of the nominal beam interaction point, 
a total reconstructed cm. energy greater than 0.5 W (yV is the T(4S') cm. energy), a 
net reconstructed cm. momentum with a z component less than 0.3 ly/c, a total energy 
deposited to ECL between 0.025 and 0.9 W ^ and a ratio Ri of the second and zeroth Fox- 
Wolfram moments less than 0.7. 

B. Dilepton Event Selection 

Lepton candidates are selected from the charged tracks by requiring Le > 0.8 for electrons 
or Cn > 0.9 and a reduced of the transverse deviation in the KLM of less than 3.5 for 
muons. In both cases a distance of closest approach to the run-dependent interaction point 
less than 0.05 cm radially and 2.0 cm in z is required. At least one SVD hit per track in 
the r-(f) view and two SVD hits in the r-z view is required. To eliminate electrons from 
7 — >■ e~^e~ conversions, electron candidates are paired with all other oppositely charged 
tracks and the invariant mass (assuming the electron mass hypothesis) Me+e- is calculated. 
If Me+e- < 100 MeV/c^, the electron candidate is rejected. If a hadronic event contains 
more than two lepton candidates, the two with the highest cm. momenta are used. 

The two lepton candidates must satisfy additional criteria. The cm. momentum of each 
lepton is required to be in the range 1.2 GeV/c < p* < 2.3 GeV/c. The lower cut reduces 
contributions from secondary charm decay. The upper cut reduces continuum contributions. 
Each lepton track must satisfy 30° < 6'iab < 135°. This cut selects tracks with better z vertex 
resolution and better lepton identification. Events that contain one or more J/ ip candidates 
are rejected. The invariant mass of each candidate lepton paired with each oppositely 
charged track (assuming the correct lepton mass hypothesis) is calculated. If the invariant 
mass falls into the J/ijj region, defined as —0.15 GeV/c^ < (Me+e- — Mj/^) < 0.05 GeV/c^ 
or -0.05 GeV/c^ < (M^+^- - Mj/^) < 0.05 GeV/c^, the candidate event is rejected. The 
looser lower cut for the electron pair invariant mass is used to reject J/tp decays with a low 
invariant mass due to bremsstrahlung of the daughter electron(s). 

As can be seen in Fig. 01 distributions of the opening angle of the two tracks in the cm. 
frame, cos 9}^, for the fifi and e/x pairs show distinct peaks in the back-to-back direction 
(cos 6'^^ ~ —1). This background is caused by jet-like continuum events and events with a 
primary lepton and a secondary lepton originating from the same B meson. Also spikes 
can be seen at cos 6'|^ = 1 in the fifi pairs. This structure is caused by jet-like continuum 
events where a non-muon track is identified as muon because it is assigned hits in the KLM 
from a true muon. The opening dilepton angle in the cm. frame 6^^ is required to satisfy 
—0.80 < cos 6*1^ < 0.95 in order to reduce this background. 

With these selection criteria there are 46551 positive and 45507 negative same sign dilep- 
ton events found in the on-resonance data. Continuum contributions are estimated to 
be 2229.8 for positive and 1556.5 for negative same sign events, based on the yield from 
off-resonance data. To estimate the continuum contribution from off-resonance data, off- 
resonance yields were scaled by the integrated luminosities and cross-section ratio. The 
scaling factor is defined in the Eq. El These dilepton yields decomposed into different lepton 
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FIG. 4: cos distributions for the dilepton samples of the on-resonance (open histogram) and 
scaled off-resonance (filled histogram) data. (a),(b) show ee events, (c),(d) show ^^u events and 
(e),(f) are from e/x combinations. (a),(c) and (e) are the ++ charge case and (b), (d) and (f) are 
the charge case. The arrows indicate the selected range. 

categories, are given in Table B 
C. Az Determination 



The z-coordinate of each B meson decay vertex is the production point of the daugh- 
ter lepton, which is determined from the intersection of the lepton track with the run- 
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TABLE I: Number of dilepton events 



combination 


on-resonance 


off-resonance 


continuum 


positive 


negative 


positive 


negative 


positive 


negative 


ee 


9059± 95.2 


9028± 95.0 


11± 3.3 


11± 3.3 


96. 2± 28.9 


96. 2± 28.9 




14672±121.1 


14014±118.4 


144±12.0 


lOOiblO.O 


1259.2±104.9 


874.4± 87.4 


efj, 


22802±151.0 


22453±149.8 


lOOiblO.O 


69± 8.3 


874.4± 87.4 


603.4± 72.6 


total 


46533±215.7 


45477±213.3 


255±16.0 


180±13.4 


2229.8±139.6 


1574.0±117.3 



dependent profile of the interaction point. |A2;| = \z{ii) — z{i2) \ is the distance between the 
2;-coordinates of the two leptons. 

In order to estimate the detector resolution in the Az determination J/tp decays to e^e~ 
and fi~^fi~ are used. In these events the two tracks originate from the same point, so the 
measured Az, after the background contribution is subtracted, yields the detector resolution. 
Candidate J/ip are selected using the same requirements as the dilepton events except the 
J/ip veto. The J/ip signal regions are defined as 3.0 GeV/c^ < M{e~^e^) < 3.14 GeV/c^ 
and 3.05 GeV/c^ < M{fi+fi-~) < 3.14 GeV/c^ and the sideband region as 3.18 GeV/c^ < 
M{i~^i^) < 3.50 GeV/c^ for both electrons and muons. 

The invariant mass distributions of J/ ip candidates are fitted to a function given by 

N{M) = hoe + hie + A{M - Bf + C. (4) 

Here h^ and hi are the heights of both Gaussians, Mq is the Gaussian mean which is common 
to two Gaussians, ctq and ai are the width of the Gaussians. A parameter S", given as S = 
for M > Mq and S = ao + a{M — Mq) for M < Mq, is introduced to modify the lower mass 
tail of one of the Gaussians for the effect of bremsstrahlung using another parameter a. 
A, B and C are the parameters of the background function. The Az distribution of the 
sideband region is scaled to the background yield in the signal region and subtracted from 
the signal region Az distribution. 

The J/ ip mass distributions and the Az distributions are shown in Fig. The RMS of 
the Az distributions are 193 (xm for J/ip ^ e'^e~ , 177 (xm for J/ip ^ A^'*'/^"; and 185 |J.m for 
the combined J/tp ^ . 



D. Subtraction of Continuum Events 

A sample of dilepton events originating from BB events is obtained by subtracting the 
luminosity and cross-section scaled off-resonance data from the on-resonance data. Since 
the kinematics of dilepton candidates are generally different in these two data samples, the 
subtraction should, in principle be performed in a six dimensional (p*, jf^, Q\, Q^, Az) 
space for each lepton flavour and charge combination, where ^^2) is the polar angle with 
respect to the beam axis of more(less) energetic lepton in cm. frame. The number of BB 
events is obtained from 

^BB^vl.vl. 0*2, o;„ Az) = NUvl.Pl ei, e*, Az) - J-^^J^n,^{pI,p;, ei, e*, 9i„ Az) 

J L-oSat Son 

(5) 
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FIG. 5: Mass distributions for J/tp — > e'^e^ (a) and J/^p — > /^^/x^ (b). The arrows indicate the 
signal region for each decay mode. The dashed hnes indicate the fitted background component and 
the soHd hnes show the total fit results. The cross hatched area shows the estimated combinatorial 
background in the signal region and single hatched area is the sideband region used to estimate 
background Az distribution. The Az distributions for J/'if) ^ e~^e~ (c) and J/^/J — > /i"*"^" (d). 
Open histograms are for all J/tp candidates in the signal region and hatched histograms are for the 
background. 

where N^g and A^on{ofr) are the dilepton yields of BB origin and on(off)-resonance data, 
respectively, / Con{oS)dt and Son(ofr) are the integrated luminosities and the square of cm. 
energies for on(off )-resonance, respectively. 6^^ is included here because this variable behaves 
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distinctly differently in the two data samples for the cases containing muons, as shown in 
Fig.ia 

Given the available statistics, this approach is not possible. Instead, we perform the sub- 
traction by weighting the on-resonance and off-resonance yields for one of the six kinematical 
variables, while integrating over the five other variables. We obtain weighting factors for the 
six variables by repeating this procedure. Since, t oa first approximation, the six variables 
are not correlated with each other, this approach provides the BB yield in the six variable 
space. The weighting factors are given by w{k) = (l/rBB){,No^{k) — f Nos{k)) / Non{k) where 
k denotes each of six variables, / is the scaling factor for the luminosity and cm. energy 
introduced in Eq. El and tbb = N^^^^/N^"^^^ is the fraction of total BB in the on-resonance 
yield after integrtating over all six variables and is used for the normalization. While the 
weighting factors show very little dependence on pi, 61, and 62 for all combinations of 
lepton flavours and charges, a clear dependence is observed for 9^^ in the case of fifi and e/i 
data samples as shown in Fig. El A clear dependence on AZ is also seen for all lepton pair 
combinations. 

Using this method, the dilepton candidate yield for each lepton flavour and charge com- 
bination is then given in terms of the on-resonance yield and the weighting factors by 

NBBipl,P*2, 01 9*, 9*u, Az) = rBBllwik)N,^{pl,pl, 9*,, 9^, 9*,,, Az). (6) 

k 

The Az dependence of the dilepton yields are obtained by projecting 
Nbb(,Pi, P2^ 01, 9^, 9}f^, Az) onto the |A2;| axis. 

IV. RESULT 

A. Corrections to lepton candidates 

The number of detected leptons for each lepton flavour and charge A^^^. is related to the 
number of true leptons by 

Ni,ip\ 01ab) = N^{P\ ^lab)4k(P*, OMe%ip*, 01ab) + J] r± ^lab)r/^,(/, ^lab)}, (7) 

h=TT,K,p 

where e^j^ and e^^^ are the efficiencies for track flnding and lepton identification, r^^ is the 
relative multiplicity of hadron h with respect to lepton i in the BB event, and r^^^ is the 
rate of hadrons h faking leptons i. The relative multiplicities are determined from B^B^ 
Monte Carlo (MC) events, and are shown in Fig. [71 

Using the measured efficiencies and fake rates and the MC determined relative multiplic- 
ities, the correction factors Ni/N^^^t were determined in 7 bins of p* /GeV/c (1.2-1.3, 1.3-1.4, 
1.4-1.5, 1.5-1.6, 1.6-1.8, 1.8-2.0 and 2.0-2.3), and 8 bins of ^i^b (30°-37°, 37°-50°, 50°-77°, 
77°-82°, 82°-lll°, 111°-119°, 119°-128° and 128°-135°). The fake rates are measured in 
the laboratory frame. For the correction, they are converted into quantities in {p*, ^lab)- 

After the correction, the dilepton sample contains true leptons that come either from 
prompt neutral B meson decay (signal) or from background processes such as charged B 
meson decay, secondary charm decay, or other leptonic B meson processes. 
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FIG. 6: cos 9^^ dependence of the weighting factor for the fraction of the dilepton yield of BB 
origin in the on-resonance data for e"*"e"^ (a) and e~e~ (b) and corresponding for ^/x (c) and (d) 
and efj, (e) and (f). 



B. Fit to Az Distribution 

A binned maximum likelihood fit with signal and background contributions is used to 
extract y4s;(|Az|) from the Az distribution. The Az distribution for the signal events is 
given by Eq. |21 assuming AF is negligible as, 

P^^ oc e^l'^^l/'so (1 _ cos(AmrfAt)), (8) 
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FIG. 7: Relative multiplicities of hadrons as a function of piab with respect to leptons. The hadron 
and lepton species are indicated on the vertical axis of each graph. Filled circles are for positive 
tracks and open circles are for negative tracks. Though in (a), (b), (d) and (e), the difference 
between positive and negative is less than 1%, in both (c) and (f), proton rate is larger than 
anti-proton rate by about 8%. 

convolved with the detector response function described earlier. Here r^o is the 5° lifetime 
and Am^ is the mixing parameter. These parameters are fixed to their world average 
values 

The backgrounds are placed into two categories; correctly tagged (CT), and wrongly 
tagged (WT). The CT category mainly contains events in which both leptons come from 
secondary charm decay in B^B^ — > B^B^{B^B^) (mixed) processes. The WT category 
contains events in which one lepton is from secondary charm decay of unmixed B^B^ or 
B^B~ and the other is from a semi-leptonic B decay. Though background Az distributions 
are estimated using MC simulations, this MC background distribution overestimates the 
Az resolution in the data. To correct for this the MC Az distribution is convolved with a 
a = (68 ±19) \im Gaussian The Az distributions for the true same-sign dilepton events 

where positive(++) and negative( ) samples are combined, are shown in Fig. |H1 together 

with the fit results. The x^/n.d.f. of the fit is 72.17/38. In the fit, the ratio of CT to WT is 
fixed at the MC value, and only the ratio of signal and background is allowed to float. The 
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MC estimated CT and WT contributions to the Az distribution are shown in Fig. |H| 



14000 




lAzl(cin) 



FIG. 8: Az distribution for the true dilepton events(++ and are combined). Points with error 

bars are data. The dot-dashed hne shows the contribution from CT backgrounds, the dotted hne 
shows the WT background contributions, the dashed line indicates the signal component and the 
solid line indicates the total of the fit. 
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C. Charge Asymmetry 



The measured same-sign dilepton charge asymmetry is defined as 

where N^^{Az) are the Az distributions of the true dilepton yields. 

Since N'^'^{Az) are the sum of signal and background, A^^^(Az) = Nf'^(Az)+N^^{Az), 
the dilepton charge asymmetry Agi is related to Au by 

, , N++{Az) - N--{Az) N,{Az) ^ , ^ , 

AuiAz) = ' ^ ... . ' — , J = Asi(Az)d(Az) 10 

where A^^ = A/'++ + N'- and = Nfj^^ + iV^"". A dilution factor, d{Az) = 
Ns{Az) / {Ns{Az) + Nb{Az)), is calculated using the signal and background yields, which 
are determined in the fit given in Fig. |S| The result of Asi{Az) which is determined from 
measured Af£{Az) and dilution factor d{Az) is shown in Fig. IHl 

The dilepton charge asymmetry is A^^dA^;!) is a time integrated quantity and does not 
depend on Az. Fitting this distribution to a constant in the region of 0.015 cm < \Az\ < 
0.200 cm, yields A^i = (-0.13 ± 0.60)% and the x^/n.d.f. is 68.73/36. The optimum fitting 
range is determined using a MC study. 



D. Cross Checks 

As a consistency check, Agi is obtained for the ee, /i/i, and efi data samples, separately. 
The results, Asi{ee) = (-1.41 ± 1.13)%, A,i{fifi) = (+1.66±1.30)%, and Asi{efi) = (-0.44± 
0.94)%, are consistent with the primary result. Here the errors are statistical only. 

In the extraction of Agi from Au{Az) using Eq. [TU[ it is assumed that N^~^ = N^^ . 
The validity of this assumption is confirmed by repeating the fit without it. This yields 
N++ = 20452 ± 126 and N^' = 20028 ± 124 in the range 0.015 cm < \Az\ < 0.200 cm, 
which is consistent with the initial assumption. 



E. Systematic errors 

Systematic errors in the determination of Agi come from uncertainties in: i) the event 
selection criteria, ii) corrections for efficiencies of track finding and lepton identification and 
for lepton misidentification, iii) the continuum subtraction, iv) the Az fit for the dilepton 
sample, v) the Az fit for the determination of Agi. 

Uncertainties in the event selection are estimated by repeating the analysis with varied 
cut values. For the track selection, the ^lab cut is varied from the nominal 30°-135° to 
51°-117° (barrel detector part only) in several steps, the closest-approach cut in r0 from its 
nominal value of 0.05 cm to 0.02 cm in six steps, the closest-approach cut in z from 2 cm to 
1 cm in five steps, the p* cut from its nominal value of 1.2 GeV/c to 1.3 GeV/c in four steps, 
the cos 9}^ cut from nominal -0.80 < cos^|^ < 0.95 to -0.70 < cos^|^ < 0.90 in several 
steps, the requirement on the number of SVD hits by ±1 from nominal value of greater or 
equal to two for z and one for rcj). In addition, the mass cuts that reject J/ip and 7 —>■ e~^e~ 
are widened by 20% and 50%. 
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FIG. 9: \Az\ distribution for Agi. 



To estimate the systematic error from the continuum subtraction, the analysis is repeated 
with the continuum subtraction varied by the statistical error of the off-resonance yield. 

The contributions from track corrections are estimated by varying each of the efficiencies, 
fake rates, and relative multiplicities by ±la. 

The contribution from the detector Az response function is estimated by changing the 
response function width according to the statistics of Az distribution of J/^ £~^£~ sample. 

The contribution from the 68 |^m smearing is estimated by repeating the analysis with 
50 [ira and 87 |xm smearing, which are the values obtained when the for the Az fit is 
changed by one compared with the default fit. The contributions from uncertainties in Am^ 
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and tbo are also estimated by varying the nominal values by ±lcr. The dilution factor fitting 
range is varied from nominal 0.00 cm < \Az\ < 0.20 cm to 0.00 cm < \Az\ < 0.05 cm. 

For the fitting range for the determination of the final Agi, the lower limit is varied from 
its nominal value of 0.015 cm to 0.040 cm in several steps. 

The results of the systematic error determination for Agi are summarized in Table IHl 



TABLE II: Source of systematic errors for the measurement of Agi 



category 


source 


AAsi (%) 


event selection 


track selection 


±0.236 




cos cut 


±0.107 




lepton pair veto 


±0.167 


continuum subtraction 




±0.314 


track corrections 


track finding efficiency 


±0.074 




electron identification efficiency 


±0.058 




muon identification efficiency 


±0.208 




fake electrons 


±0.031 




fake muons 


±0.047 




relative multiplicity 


±0.057 


Az fit for dileptons 


detector response function 


±0.009 




Arud 


±0.011 




Tbo 


±0.009 




68 (J.m smearing of background Az 


±0.009 




statistics of background MC 


±0.022 




fitting range 


±0.004 




assuming -/V^^ = 


±0.139 


Az fit for Asi 


fitting range 


±0.207 


total 




±0.561 



V. CONCLUSION 

The charge asymmetry for same-sign dilepton events T(4S') decays has been measured. 
The result is related to a CP- violation parameter in B^-B^ mixing, Agi = (—0.13 ± 
0.60(stat) ± 0.56(sys))%, or equivalently \q/p\ = 1.0006 ± 0.0030(stat) ± 0.0028(sys). The 
measured Agi is consistent with zero, or equivalently, \q/p\ is consistent with unity. This 
implies CP- violation in B^-B^ mixing is below the 0(10-^) level. The CP-violation param- 
eter es can be calculated as Pe(eB)/(l ± |eBp) = (—0.3 ± 1.5(stat) ± 1.4(sys)) x 10~^, using 
the exact formula 



l±v/(l-^/)/(l±^0' 

These results are consistent with and provide significantly more restrictive bounds than 
previous measurements [11 
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